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Abstract 

Electrochemical Scanning Tunneling Microscopy (ECSTM) is gaining popularity 
as a tool to implement proof-of-concept single (bio) molecular transistor. The under- 
standing of such systems requires a discussion of the mechanism of the electrochemical 
current gating, which is intimately related to the electrostatic potential distribution 
in the tip-substrate gap where the redox active adsorbate is placed. In this article, we 
derive a relation that connects the local standard potential of the redox molecule in 
the tunneling junction with the applied electrode potentials, and we compare it with 
previously proposed relations. In particular, we show that a linear dependence of the 
local standard potential on the applied bias does not necessarily imply a monotonous 
potential drop between the electrodes. In addition, we calculate the electrostatic po- 
tential distribution and the parameters entering the derived relation for ECSTM on 
a redox metalloprotein (Azurin from P. Aeruginosa) , for which experimental results 
exist. Finally, we give an estimate of the gating efficiency when the ECSTM setup 
including Azurin is interpreted as a single biomolecular wet transistor, confirming 
the effectiveness of the electrochemical gating for this system. 



1 Introduction 



The Holy-Grail of molecular electronics^ is the realization of working transistors based on a 
single molecule. In particular, to realize a field effect single (bio)-molecule based transistor, 
three different electrodes should surround the molecule: a source, a drain and a gate. Ide- 
ally, source and drain should exchange electrons with the molecule (i.e., they should provide 
the electric current) without affecting the molecular energy levels, while gate should tune 
the molecular energy level (via the produced electric field) without exchanging electrons 
with it. This implies that, when the gate is implemented as a third metallic electrode, it 
must be microscopically close to the molecule and, at the same time, electrically isolated 
from it.^ Meeting such requirements represents a formidable experimental task. An el- 
egant alternative solution to the gating problem is provided by the electrochemical gate 
effect,^ conveniently implemented with Electrochemical Scanning Tunneling Microscopy 
(ECSTM).^ 

When standard (i.e., non electrochemical) STM is used to study single redox-active 
(bio)molecule adsorbed on a solid conductive substrate, the resulting setup is close to that 
of a field effect transistor (the conductive substrate represents the source and the tip acts 
as the drain) but the gate electrode is lacking. In ECSTM, the STM tip and substrate 
are immersed in an ionic solution, and their potentials are controlled with a bipotentiostat 
w.r.t a reference electrode. Thanks to this bipotentiostat, ECSTM allows to control not 
only the potential bias between tip and substrate, as in in vacuo STM, but also another 
parameter: the difference between the substrate (tip) potential and the reference potential 
of the solution. In other words, in ECSTM we can modify the electrochemical potentials of 
tip and substrate with respect to that of a redox molecule in solution, or, vice- versa, we can 
tune the energy of the molecular redox level w.r.t. the tip and substrate electrochemical 
potentials. Since the relative position of the redox level controls the molecular conductance. 
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gating is eventually obtained. The features of ECSTM has been exploited by different 
groups to demonstrate electrochemical gating,^"^^ and thus to propose/implement proof- 
of-concept single molecule transistors.^^' 

Several theoretical works have been devoted to the study of the conduction mechanisms 
of redox adsorbate in ECSTM^^"^° and, in general, to STM in water solution. ^^"^"^ In the 
present article, the focus is not on the conduction mechanism itself, but on the electro- 
static potential induced by the applied electrode potentials in the tip-substrate gap and 
on its relations with the electrochemical gate effects. A fundamental role in interpreting 
experiments on standard STM and molecular wires is played by the electrostatic poten- 
tial distribution in the tunnehng gap. It has been argued that the feature of the I{V) 
curve dramatically depends on the fraction of tip-substrate potential bias dropping on the 
molecule.^^ In ECSTM, it is clear that electrostatic potential must play a relevant role too. 
This was recognized early, however, till very recently,^^ not much attention has been given 
to the effects of the ionic atmospheres on the electrostatic potential distribution induced 
in the tip-substrate gap by the applied electrode potentials. For example, in interpreting 
experiments on adsorbate conduction in ECSTM, a hnear relation between the fraction 
of the bias potential acting on a redox center and the geometrical position of the redox 
center in the tunneling gap has been assumed. Is this picture correct for ECSTM of 
adsorbates? One of the aims of the present article is to answer to this question in a specific 
case that is particularly ambiguous. 

In addition, the feature of the electrostatic potential distribution in ECSTM are strictly 
related to a quantity which is fundamental when considering ECSTM experiments as re- 
alizations of single molecule transistors: the gating efficiency. In this article, this quantity 
will be estimated as a function of the system parameters. 

The discussion of electrochemical phenomena such as those controUing gating in EC- 
STM may become muddled due to the use of the term potential to indicate different quan- 
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titles. In the discussion that follows we should be careful in distinguishing between these 
quantities: the electrostatic potential (physical dimension of an energy divided by a 
charge); the electrochemical potential fi (physical dimension of an energy) and the electrode 
potential EoiV (physical dimension as 0, but it is related to a difference of electrochemical 
potentials, divided by e, the modulus of the electron charge). 

The rest of the article is organized as follows: in Sect 12] we review the basics of the 
gating mechanism in ECSTM, beginning with the ideal case and describing the necessary 
refinements associated with the microscopic origin of the process; in SectJS] we present 
the continuum model that we use to describe the potential in the ECSTM setup while 
in Sect JH results are presented and discussed. In particular, a simple relation between the 
standard potential of the redox molecule in the gap and the applied tip/substrate potentials 
is derived, the electrostatic potential distribution for a protein in the tip-substrate gap^ is 
discussed and the efficiency of an ECSTM transistor based on Azurin is presented. Finally, 
in Sect J5] some conclusions are drawn. 

2 Review of the gating mechanism 

In this section, we present a short review of the electrostatic aspects of the gating mech- 
anism. A description of the gating mechanism in the framework of ECSTM electron tun- 
neling mechanisms can be found, e.g., in Ref.28. 

We start our discussion by considering an ideal case. Let us suppose that the substrate 
(i.e., the source in the transistor nomenclature) is grounded. Thus, all the electrochemical 
potentials /i are naturally referred to that of the substrate, figub- In other words, imposing 
a variation to the substrate-reference potential difference is equivalent to change the elec- 
trochemical potential of the reference electrode (/i^e/)-^^ When such a variation is imposed 
(i.e., when firef is changed), the characteristic semi-reaction of the reference electrode is 
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drawn out of thermodynamic equilibrium. To recover the equihbrium, the reaction ad- 
vances or retrocedes, releasing or removing ions from the solution. This change in the ionic 
balance of the solution modifies (psou the electrostatic potential of the solution itself, and, 
in turn, the electrochemical potential of the ions already present in solution. The new 
equilibrium is settled when the electrostatic potential variation of the solution {^(psoi) is 
such that to have compensated for the initially imposed electrochemical potential varia- 
tions: — eA0so« = ^fJ^ref, where e is the modulus of the electron charge (see appendix A for 
a proof ).^° The molecule feels now the new (psoi- The energy of each molecular electronic 
level (in particular of the redox one, i.e., that occupied upon reduction) is thus modified by 
a quantity —e^(f)soi = AyU^e/- In conclusion, if the potential difference applied between the 
transistor source (i.e., the ECSTM substrate) and the gate (i.e., the reference electrode) 
is changed by a given quantity, then, in this ideal picture, the energy of the redox level of 
the molecule is changed by the same amount (times — e), while the source (substrate) and 
drain (tip) remained constant. This modification of the redox energy level of the molecule 
w.r.t. the source (substrate) and drain (tip) level is precisely what gives origin to gating. 

So far, we did not stress that while (psoi varies in response to a change in firef, the 
electrostatic potential in the substrate is fixed, since such electrode is grounded. Thus, 
the electrostatic potential difference between the interior of the substrate and the solution 
must change of A0so« = A/ij.e//(— e). Ideally, there is a step in the potential profile going 
from the substrate to the solution, whose size is equal to the applied substrate-reference 
potential up to an unknown but fixed constant. This electrostatic potential drop is step- 
like only at the macroscopic scale (see Fig. [1]). In fact, microscopically, the electrostatic 
potential drop is due to the imbalance of positive and negative ions close to the charged 
electrode surface, which creates a net density of charge in the solution layer neighboring 
the electrode. The structure of this density of charge is described by the Stern double- 
layer model, which extends for a few nm's (depending mainly on the ionic strength) in 
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Figure 1: Schematic representation of the electrostatic potential drop from the substrate 
to the solution assuming an ideal (left) or real (right) behavior. Note that when (psoi is 
changed, the electrostatic potential at the position of the redox center changes by different 
amounts in the two cases. 

the solution. When the redox molecule resides inside this layer (as for adsorbates) the 
electrostatic potential acting on the molecule (called 0^) is different from (psoi- Thus, the 
electrochemical potential variation imposed on the gate (i.e., reference) electrode does not 
translate completely into the variation of (pm- A measure the effectiveness of the gate in 
controlling the molecule potential is the gate control parameter,"^ defined as: 



P 



d(j)„ 



(1) 



where Vg is the potential applied to the gate. Clearly, /5 = 1 in the ideal case while [3 <1 
in the real case. 

We end this discussion by remarking that the effects related to the incomplete screening 
of the electrode charge at the position of the molecule have well-known consequences in 
electrochemical kinetics, where they originate the so-called Frumkin effect. ^^'^^ 

In the following, we shall discuss how the electrostatic effects of the source and drain 
electrodes on (pm (and thus on (3) can be calculated within a simple model. 
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3 Outline of the model 



We are not aware of calculation of electrostatic potential distribution for a molecule (in par- 
ticular, a protein) in ECSTM, although calculations for empty tunneling gap appeared very 
recently.^^ Thus, we start by considering a relatively simple model based on a continuum 
dielectric representation, for which useful relations can be promptly derived. Basically, 
tip and substrate arc represented as perfect conductors (substrate with a planar surface, 
tip with a spherically curved surface), and the adsorbate is a dielectric with a realistic 
shape (solvent-excluded surface) . The redox center of the adsorbate is not exphcitly mod- 
eled: we simply define (pm as the electrostatic potential at the geometrical position of such 
center. As for the ionic solution surrounding the adsorbate, it is described as another di- 
electric while the ionic atmosphere is implicitly accounted for by considering the linearized 
Poisson-Boltzmann equation to obtain the potential in the system: 

_V-(eV0)+eA;2(0-0,,,) = ^ (2) 

where e is the relative dielectric constant (of the solvent or the adsorbate, depending on 
the position) and k is the inverse screening length (equal to zero inside the adsorbate), 
proportional to the square root of the ionic strength I. The Helmholtz layer is taken into 
account by assuming that the region of solution closest to the electrodes cannot contain 
ions (i.e.. A; = 0, see appendix B). The electrostatic potential drop inside this region is 
thus linear. We shall not consider any specific ion adsorption at the electrode surfaces, 
i.e., the concentration of ions in the Helmholtz layer is only related to their electrostatic 
preferences and not to other favorable interactions with the surface. We remark that a 
reduced dielectric constant for the water layer in the close proximity of the electrodes^^ 
has been used as discussed in appendix B. 

It is clear that the present model is quite simplified and neglects a number of effects, 
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such as the non-hnear character of the Poisson-Boltzmann equation, specific ion adsorption 
and the accounting of the ion finite size in ion-ion interactions. Nevertheless, it catches 
the main physical features of the ion screening effects (also in complex systems involving 
biomolecules) and, more important than this, represents the framework in which corrections 
to the molecular redox level energy as simple as that used in Ref.l9 can be derived. Refined 
calculations based on the non-linear Poisson-Boltzmann equation are given and discussed 
in Sect. 14.31 for a specific system. 

The electrostatic potential that we consider is only due to the applied potential dif- 
ference: the contributions from the intrinsic charge distribution of the redox adsorbate 
is already included in the factors that determine the standard redox potential; charging 
effects of the molecule due to the electronic coupling with the electrodes are not consid- 
ered here, since we focus on the electrostatic effects due to the charged electrodes. In any 
case, in the framework of our simple, linear model, these effects are superimposable to our 
results. 

4 Results and Discussion 

4.1 The relation between the electrostatic potential acting on 
the molecule and the potentials applied to the tip and the 
substrate 

In this section, we present a general relation between the electrostatic potential acting 
on the molecule redox center (0^) and the potentials imposed on the electrodes {Eup 
and Esub)- Such a relation is based on the continuum model presented above. First 
of all, we have to specify the electrostatic boundary conditions at the tip and substrate 
surfaces. We assume that the electrostatic potential at the substrate surface (psub is given 
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by (psub = Esub — E^ub + '^soi, where Egub is the potential imposed to the substrate w.r.t. 
the reference electrode while E^^l is the potential of zero charge (pzc), i.e., the substrate 
potential at which the surface charge is nil. Analogous definitions are used for the tip, 
yielding (p^p = Eup — E^^^ + (psoi- We remark that the values of E^^l and E^^^ refer to a 
system where tip and substrate are well separated. Clearly, this does not mean that, e.g., 
the substrate involved in a STM has a null surface charge at E^^, since it is exposed to 
the electrostatic effects of the other electrode. This is naturally taken into account by the 
present model. 

Since the equation that defines the potential in our model, eq.([2]), is linear in — (^gou 
one obtains: 

4>m — 4>sol = atip{,4>tiv ~ (f'sol) + Oisubi4'sub — 4>sol) = 0!tip{Etip — E^^^^) + OsubiEsub — -^^uD (3) 

Thus, the effect of the tip and the substrate potentials on the redox center is, in our 
simple model, completely described by an intrinsic electrode property, E^^'^, and by the 
two coefficients aup and agub- When there is no electrostatic coupling between tip/substrate 
and the redox center, aup = asub = 0; when the redox level is pinned to, e.g., the substrate 
then asub = 1- Clearly, < atip,asub < 1- In addition, it is possible to proof that 
C(tip + C(sub < 1- Before giving the results of the actual calculations of asub and aup, we 
would like to draw a connection with another relation described in the literature to take 
into account the effects of the tip and substrate potential on the molecular redox levels. To 
account for such effects in their most recent models, ^^'^^ Ulstrup, Kuznetsov and coworkers 
substitute the overpotential Esub — E^ (where E'^ is the standard redox potential for the 
molecule) with an effective overpotential given by C,{Esub — E'^) — •yVuas with < 7 < 1. 
In ref.l9b, ^ and 7 are defined as "the fractions of the overpotential {^) and bias voltage 
(7) at the redox site" . Following our arguments, in the ECSTM environment the effective 



9 



redox potential E^ i^^ is given by i?^ + 4>m — 4>soh thus the effective overpotential becomes 
Esub — (-f'm + 0m — (psol)- Equating this expression with that used in Ref.l9 leads to: 

Esub - {El + <f>m- <Psol) = CiEsub - El) - -iVuas (4) 

which means: 

0m - (Psol = (1 - - El) + -iVuas (5) 

In order to compare eq.(l5]) with eq.(l3]), we can recast the latter in the form: 

0m — (Psol = {C^tip + Otsub) {Esub " -E'm) + CltipVuas + + «snfe) E^ — UupE^^^ — UsubElH (6) 

By comparing eq.© and eq.([6]), one can identify ^ = 1 — {aup + Osnfe) and 7 = aup. 
However, there are different terms in eq.(l6]) that does not appear in eq.(l5]), playing the 
role of a correction to the molecular standard redox potential E^. Not surprisingly, such 
a correction depends on E'^^'^ of the tip and the substrate. Notably when ,^ = 1,7 must be 
zero since ^ = 1 implies aup = asub = 0. 

We remark that 7 = aup is not directly related to the drop of the applied bias in the 
tunneling gap, and thus does not reflect the position of the redox center. In fact, 7 = aup 
may be small even when the molecule is closer to the tip than to the substrate, if the 
ionic concentration is high enough (see below for an example). This is simply understood 
when one consider the peculiarity of ECSTM w.r.t gas-phase STM. In fact, when STM is 
performed in gas-phase, the electrostatic potential due to the applied bias is often (and 
reasonably) assumed to linearly (or at least monotonically) change between the tip and the 
substrate. In ECSTM, the situation is very different, as demonstrated in Ref.27. Suppose 
that tip and substrate are far enough to allow a complete screening of their respective 



10 



(b) 



tip ^ sub 

Figure 2: Electrostatic potential profiles in pure ionic solution for two different combina- 
tions of applied tip and substrate potentials, (a) Egub > E^ub^ ^up > E^^^] (b) Eg^h > E^^l, 
Etip < E^^p. Parameters of the calculations: Ionic strength I = 0.05 M, tip-substrate 
distance: 40 A 

surface charges by the ionic atmosphere in solution. Depending on these charges (which 
depend, in turn, on the imposed electrochemical potentials), we can have four different 
situations, two of which are sketched in Figj2] as examples. It is evident that the existence 
of a monotonous behavior depends on the applied potentials. What is important to stress 
here is that even when the potential drop is not linear, a linear relation between E^i^^ 
(the standard potential of the molecule in the precise ECSTM gap) and the applied Es^f, 
and Vbias is obtained from eq.(l6]): 

E'm,ioc = El + {aup + a..6) {E,^b - E^) + aupVuas + E^"^^ (7) 

This means that care must be taken to assume a linear potential drop in the interpretation 
of ECSTM experiments only on the basis of linearity of the experimental results with the 
applied potentials. 

In eq.(I7j), E''"""^ is a correction term that does not depend on Esub and Vuas- It depends 
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on E^lp and (whose experimental estimate may be difficult) but it can also collect 
variations of -EmZoc ^^^^ does not originate from the electrostatic effects that we are ex- 
ploring in this article (due, e.g., to changes in the adsorbate related to the interaction with 
the electrodes or with other adsorbates or to the different properties of water confined in 
the nanometric tip-substrate gap). 

-^m loc is important quantity in the theories of adsorbate conduction in ECSTM, since 
it is related to the position of the peak of the current (or apparent height) as a function 
of the gate potential. Thus, eq.([7]) is not only useful for theoretical considerations, but 
also as a tool to interpret the experimental data. Remarkably, when the ion screening is 
ineffective (i.e., much larger than the gap size), one obtains aup + asub ~ 1 (^ = 0) 
and - El^,^^ = -atipVuas - E'^'^^'^ (or Eup - E^^^^ = (1 - aup)Vuas - E^"^'), i.e., the 
relations between the potentials of the electrodes and the effective energy of the molecular 
electronic levels becomes similar to that used for gas-phase (or non- ionic) STM.^^ 

4.2 Electrostatic potential distribution in ECSTM experiments 
on Azurin 

When the solute is relatively large, like a protein, the electrostatic potential distribution 
in the tip-substrate gap is more complex than what depicted in Fig 121 In fact, the space 
occupied by the solute is not accessible to the ions, which cannot thus screen the potential 
due to the surface charges of the tip and the substrate. Thus, the potential drop inside 
the protein will be intermediate between two limit cases: (a) equal to that in the pure 
ionic solution; (b) similar to that in a non- ionic dielectric, i.e., linear as in gas-phase. One 
may thus wonder whether in real experiments, the potential drop is more similar to that 
in gas-phase or that in ionic solution without the adsorbate. To address this question, we 
have calculated the potential distribution with the numerical model described in appendix 
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Figure 3: Pictorial representation of the Azurin orientation in the ECSTM setup exploited 
in Ref.Gb. 

B. In particular, we choose to focus on the recent experiments by Facci and coworkers.^'' 
They studied via ECSTM the transport properties of tip-immobilized Azurin, a blue-copper 
protein, and they demonstrate a transistor like behavior for the system. In this experiment, 
Azurin was anchored to the gold tip by exploiting the accessible disulfide bridge between 
Cys3-Cys26. In our calculation, Azurin has been oriented with the major axis of inertia 
perpendicular to the substrate, with the S atoms of Cys3 and Cys26 at bonding distance 
from the tip surface. The protein redox active site (which comprises a redox active copper 
ion) is on the opposite side of the protein. A pictorial representation is given in Fig|3l This 
orientation is compatible with the experimental morphology data of Alessandrini et al.^^ 
In FigHlwe report the potential distribution along the axis perpendicular to the electrodes 
and crossing the protein. To give an idea of the complexity and inhomogeneity of the 
potential distribution we also report, in the same figure, the potential on a slice in the 
protein region. This figure confirms what we anticipated above: the potential distribution 
in the protein region is intermediate between the gas-phase and the pure ionic solution 
cases, and it is closer, in the present case, to that of a pure ionic solution. In particular, 
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Figure 4: Potential distribution for tip-immobilized Azurin in ECSTM. (a) Color-code 
figures of the potential on a planar surface perpendicular to the electrodes surfaces, (b) 
Electrostatic potential along the axis perpendicular to the electrodes and crossing the 
protein (dashed white line in panel (a)). For comparison, the electrostatic potential cal- 
culated in ionic solution without the protein ("ions, no protein") and for an homogeneous 
dielectric filling the space between the electrode ("no ions") are given. In (b), Azurin 
occupies the z region on the right of the vertical red line. Parameters of the calculation: 
Eu^ = -0.3 V+Ef/, Esub = -0.2 Y+EZl (thus Vbias = -0.1 V), ionic strength / = 0.05 
M, tip-substrate distance: 40 A, tip curvature radius: oo. 
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Figure 5: Position of the redox center in the tip-substrate gap [zm)-, estimated on the basis 
of the 7 = value. "Calc" refers to our electrostatic calculations, "Lin" refers to a 
simple linear interpolation. Parameters for the calculation: ionic strength / = 0.05 M, 
tip-substrate distance: 40 A, tip curvature radius: oo. 

we numerically show that the STM-like picture (linear potential drop between the tip 
and the substrate) is inappropriate to describe the potential drop in ECSTM even for a 
large ion-excluding adsorbate. These results also suggest that the correlation between the 
geometrical position of the redox center and the potential acting on it is not straightforward. 
To better circumstantiate this point, we have taken the hypothetical point of view of 
someone who measured the 7 = aa^ coefficient from an experiment (e.g., fitting the current 
as a function of the potentials E^uh and Vuas)-, and who wants to use this quantity to 
estimate the geometrical position of the redox center in the tunneling gap (named Zm)- 
If one assumes that Vuas drops linearly, the position will simply be given by Zm = 
7 • Igap, where Igap is the gap length (i.e., the tip-substrate distance). We can compare this 
prediction with what can be obtained by our electrostatic calculations (these calculations 
give 'y{zm), which can be inverted to yield Zmi'y))- The comparison is performed in Fig. [51 
As it can be seen, the differences can be substantial, and the redox site can be misplaced 
by the linear interpolation up to 16 A, which is the 40% of the tip-substrate gap size. It 
is also evident that the linear interpolation results always place the redox center closer to 
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the tip than what it really is. 

4.3 Effects of non-linearity in the Poisson-Boltzmann equation 
on the potential distribution 

As mentioned above, the relations and the results presented in Sects. 14.11 and 14.21 relay, 
among other assumptions, on the linearization of the Poisson-Boltzmann equation. Due 
to the potential range (0.2-0.3V) and ionic strength (0.05M) considered, this assumption 
needs to be check. In the next section, we will show that experimental results for Azurin are 
indeed linear in the applied potentials, indicating that non-linear effects in the real system 
should be small. Nevertheless, as a consistency check of the model, we have performed 
some calculations with the non-linear Poisson-Boltzmann equation to give an estimate of 
the importance of non-linear effects for the studied system. First, let us note that in 
principle, when the assumption of linearity is removed, we cannot consider separately the 
potential due to the intrinsic density of charge of Azurin and that due to the charged 
electrode surfaces. However, the total charge of Azurin at neutral pH is very small (0 or 
-1, depending on the oxidation state of Cu),^^ and electrostatic potential calculations shows 
that the molecule produces a potential in solution smaller than fc^T, that is well described 
by the linearized PB equation. Thus, to give an estimate of the effects of non-linearity 
we have focused on the potential generated by the charged electrodes alone. In particular, 
we have repeated the 3D calculations leading to FigH] by exploiting the non-linear PB 
equation. Results are presented in FiglHl As it can be seen, the general appearence of 
the two figures is very similar, pointing to the fact that the linearized Poisson-Boltzmann 
equation correctly catches the main features of the potential distribution. The color-coded 
potential distributions of the two panels (a) can be hardly distinguished. Note in particular 
that the relative behavior of the protein vs the ionic solution potential drop is the same in 
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Figure 6: Potential distribution for tip- immobilized Azurin in ECSTM calculated from the 
non-linear Poisson-Boltzmann equation, (a) Color-code figures of the potential on a planar 
surface perpendicular to the electrodes surfaces (the color scale is the same as in FigjH 
(b) Electrostatic potential along the axis perpendicular to the electrodes and crossing the 
protein (dashed white line in panel (a)). The electrostatic potential calculated in ionic 
solution without the protein ("no protein") is also given. 



17 



-0.3 



tip 10 20 30 sub 

z(A) 



Figure 7: Electrostatic potential profile in the tunneling gap calculated in ionic solution 
without the protein and neglecting the Helmholtz layer effects, "lin" refers to the potential 
obtained from the linearized PB equation while "non-lin" refers to the non-linear PB 
results. 

the two figures. In panels (b), the most evident difference is in the region between and 7 
A, outside the protein. On the absolute potential scales, a difference of 15-20 mV ( 5 — 10% 
of the applied potentials) between the linear and the non-linear calculations is obtained on 
the the precise values of the potential plateau reached in the middle of the gap. The overall 
agreement between FigJH and Figj6] can be explained by noting that most of the potential 
drop takes place in the Helmholtz layers that, being free from ions, behaves similarly in the 
linear and non-linear cases. Such an idea is supported by calculations that we performed 
for the empty gap by neglecting the existence of the Helmholtz layers. The results of such 
calculations are reported in FiglTJ where the linearized and the non-linearized PB potential 
drops are plotted together. As it can be seen, linearization is a much worse assumption in 
this case, leading to larger discrepancies between the potential profiles. 
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4.4 Comparison with experimental results 

Comparing the theoretical predictions reported in this article with existing experimental 
results is not an easy task. In fact, the experimental i{Esub, Hias) curve depends on other 
parameters beside agub, C(tip and E"^"^^ that appear in eq.©, such as the reorganization 
energy and the transfer integrals of the electrode-molecule electron transfer reaction. In 
principle, to extract all these quantities one can assume a theoretical relation i{Esub, Vtias) 
and then adjust the parameters entering such a relation {agub, c^up, etc.) to reproduce the 
experimental results. However, such a multidimensional fit is not straightforward (different 
sets of parameters can give similar fitting quality). Moreover, different mechanisms have 
been proposed to explain the conduction of redox-active molecules,^^'^*^'^^ yielding different 
relations and thus different fitted parameters. Nevertheless, we have exploited the exper- 
imental data of Alessandrini et al.^^ to extract one quantity that can be compared with 
our theoretical calculations. We have chosen this experiment among the few available ones 
for a number of reasons: (a) the ionic strength is not too high (0.05M), which improves 
the quality of some of the model approximations (e.g., linearization of Poisson-Boltzmann 
equation); (b) it is one of the few single-molecule ECSTM experiments in which the current 
as a function of the substrate potential is collected for several different bias potentials; (c) 
it is the system where the non-linear nature of the potential bias drop is less obvious, since 
the ion screening is ineffective inside the protein and the ionic strength is relatively low; 
(d) the electronic couplings of the electrodes with the redox center is small, thus justifying 
the neglecting of electrostatic potential changes due to the effects of the flowing current 
on the molecular density of charge; (e) the configuration of the molecule with respect to 
the electrode is relatively well characterized, since azurin has only one chemisorption site 
available for Au (the Cys3-Cys26 solvent exposed disulphide bridge), and it is a quite rigid 
protein; (f) the size of the tip-substrate gap is, very likely, comparable to the molecular 
size ( 4 nm) and thus relatively large. Hence, ions concentration in the gap should not be 
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affected by steric exclusion plienomena.^'' 

Tlie results of Ref. 6b are curves i{Eup) at different values of Hias- Each i^Eup) curve 
has a maximum for a given value of E^p, called E^^^. Correspondingly, we can define E'^^ , 
which is simply given by E^^"" — Vbias, and -E^^^ as (-E^™''' + E^'^g^ represents 

the mid-point between the potential of the tip and that of the substrate and, in different 
theories proposed to explain the conduction of adsorbate, it is directly related to E^ (see 
Ref.20 and references therein). In addition, it does not privilege any of the two electrodes, 
and we shall use it for our theory-experiment comparison. Extracting the value of E]^^^ per 
each value of Vi,ias in the experiment, we end up with an experimental E^^^ vs Vuas curve. 
On the other hand, from eq.(I7j) it is possible to derive a t/ieoretzca/ relation E^1j^{Vbias) that 
can be used to fit the experimental E^'^^^ vs Vbias trend by using a^p and asub as adjustable 
parameters. The quality of the fit will be a measure of the validity of eq.(I7]), and thus of 
the linearization of the Poisson-Boltzmann equation, while the fitted values of asub and 
atip will be compared with the estimates based on our numerical model. To derive such a 
theoretical equation, we need a relation between E'^ i^^ and E"^^^ . The sequential two-step 
model (the most likely mechanism for the experiments in Ref. 6b, see also 20) predicts 
Em,ioc ~ E'^vg-^^ Introducing such a relation in eq.([7]) leads to the E'^^^iy^as) expression 
that we were looking for: 

1 ^ ^ rpcorr 

rpmax ^ ^tip ^sub pO j ^ (q\ 

Eavg = 9 1 — ^bias + J^m + " — 

^ J- O^tip O-sub J- (J^tip O^sub 

Fitting the data of Ref.6b, we obtained that the experimental relation between -E™^^ and 
Vbias is indeed linear as predicted by eq.® (regression coefficient R=0.95) and that the 
slope {atip — asnfe)/2(l — aup — asub) is -0.12. From our electrostatic calculations, we found 
asub = 0.26 and aup = 0.13, yielding {aup - asub)/'2{l - aup - asub) = -0.11, in good 
agreement with the experiment. Unfortunately, the intercept of eq.([8]) depends on E'^°^^ , 
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which cannot be directly derived by our simple numerical model. 

Clearly, this theory-experiment comparison is quite limited and involves some assump- 
tions (such as the choice of the conduction mechanism). However, we can at least say that 
available experiments seem not to contradict our predictions. We hope that our results 
may stimulate further experiments in this field. 

Other experiments on Azurin,^^ which employed different ECSTM geometries, cannot 
be treated with the same procedure used for the data of Ref. 6 since they do not measure 
currents but only changes in apparent heights. In addition, they do not explore the Vbias 
dependence of the results. 

4.5 Quantifying the gating efficiency /3 for the Azurin-based EC- 
STM transistor 

We now pass to the discussion of the gating efficiency for Azurin in ECSTM. We have 
already mentioned the gate control parameter /3, which is a synthetic way to express the 
gating efficiency in a molecular Field Effect Transistor (FET)-like geometry. When (3 = 1, 
all the potential applied to the gate is felt by the molecule; vice-versa when (3 = the 
gate is not able to modify the redox energy level. When ECSTM is seen as a FET, we can 
identify the potential of the gate Vg in eq. ([H) with that of the reference electrode, measured 
w.r.t the grounded electrode (the source), meaning that Vg = {firef — fJ^sub) / {—e) ■ Since 
Esub is defined above as the potential difference between the substrate and the reference 
electrode, and the substrate electrochemical potential is fixed (we recall that the substrate 
is grounded), we can write E^ub = -Vg. Thus, Eup = Esub + Hias = -Vg + Vuas- In 
addition, we have noted above that —eA(f)soi = ^f^ref, which means that the potential of 
the solution 0so; differs from by a fixed but unknown constant c (i.e., ((>soi = Vg + c). 
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Substituting these relations in eq.([3]), one obtains: 



(t)m = c + Vg + aup{-Vg + Vuas ' Ef^ + a,ub{-Vg - EZl) (9) 
j3 can be easily derived from eq. Q : 



= 1 - {atip + asub) (10) 



The meaning of eq. (fTO|) is clear: when the redox center does not feel the surface charge of 
tip and substrate (i.e., aup = asub = 0), the gating works ideally (/? = 1); on the contrary, 
if the molecule is affected by the tip and/or the substrate charge the gating efficiency is 
reduced (0 < /5 < 1). Remarkably, the value of /? also affects the position of the peak 
in the tunneling current as a function of the gate potential. In fact, the peak position is 
related to the local standard potential of the adsorbate, E'^ i^^, which depends on P via 
atip and agub as illustrated by eq.©. For example, in the sequential two-step model the 
substrate potential giving the current maximum, E^^, satisfies E^^ + Vbms/2 ~ -E'mZoc-^° 
Introducing such condition in eq.([7]) one obtains: 

^ 1/9 rpcorr 

Esub -Em^ p Vbias + (iij 

When P = 1 (i.e., aup = a^ub = 0), then E^T = + - while /5 < 1 can 

either increase or decrease E^l^ depending of the sign of E'^°^^ and the value of aup. 

Since an effective screening of the electrode surface charges is needed to have high values 
of /3, one expects that an important parameters for determining f3 is the ionic strength I of 
the solution. In fact, the higher the ionic strength, the smaller is the screening length k~^. 
To estimate (at least in the framework of our simple, linearized model) the dependence of 
P on k~^, we have performed calculations on the Azurin-based biotransistor, whose results 
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Figure 8: The gating efficiency (3 for an Azurin-based transistor as a function of the 
screening length k~^. The tip-substrate distance is 40 A 

are depicted in FigJHl Although the values of (3 indicate that the azurin-based transistor is 
not working ideally, they are high enough to con&m that the ECSTM approach to gating is 
convincing even in the ambiguous situation of a protein. In particular, the values of /? range 
from 0.45 to 0.85 for k^^ between 25 A and 5A. The value of (3 for the ionic strength used 
in the experiment of Ref. 6b (1=0.05 M) is around 0.6, which is coherent with the observed 
gating behavior. In the light of the results of Fig. [HI we can define k~^, the typical length 
over which the effects of a charge in ionic solution are screened, as an effective gate size 
in ECSTM. In fact, (3^1 when k^^ << Igap and /5 ~ when k~^ » Igap, as expected. 
Remarkably, we also found that (3 = 0.5 when k~^ ^ 20A ^ ^gapf^ (the tip-substrate 
distance Igap is 40 A, which is in turn related to the protein size). 

5 Conclusion 

To summarize the contents of this article, we have discussed how the gating mechanism 
acts in a wet transistor based on ECSTM and we have used a simple computational model 
to give numerical estimates of the relevant quantities for a bio-transistor based on Azurin. 
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In particular: 

• We argued that the potential drop (due to the apphed bias) in ECSTM of large ad- 
sorbates is not monotonic and linear, in analogy to what happens for empty ECSTM 
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gaps. 

• We showed that the potentials applied to tip and substrate can indeed affect the redox 
energy level of the probed biomolecule, when realistic parameters for the system are 
used. 

• We demonstrated that, in the simplifying framework of a linearized Poisson-Boltzmann 
model, the variation of the redox energy level is indeed proportional to the tip and 
substrate potentials, as exploited in models to explain ECSTM. We have empha- 
sized that this proportionality does not generally imply a spatially linear potential 
drop between tip an substrate, an assumption that may lead to misinterpretations of 
the experimental results. In addition, we have considered the role of pzc potentials 
in a straightforward way. 

• We have quantified the gating control parameter (3 for the protein bio-transistor, an 
important quantity that expresses the gating efficiency as a function of the parameters 
of the system. 

• We have identified a quantity (the screening length k~^) that acts as an effective gate 
size in ECSTM. 

We remark again that the computational model is quite rough, and to obtain more precise 
estimates, more complex models must be used. However, since this is, to the best of our 
knowledge, the first study on the potentials in these bio-ECSTM systems, the simple model 
used allowed a clear indication of trends and gave insights into the basic physics of the 
system. 
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A Proof that —eAcpsoi = AfJ^ref 

Let us suppose that the half-reaction of the reference electrode is in the general form: 

J2u,r^ -ne;i^^ = (12) 
J 

where we have used the convention that stoichiometric coefficient are negative for reac- 
tant and positive for product, while cj is the charge of the species J and n is the number 
of exchanged electrons. Charge balance of eq. (fT2!) implies: 

Yl ^j^i + n = Q (13) 
i 

while thermodynamic equilibrium is expressed by: 

l^j^-j - rifiref = (14) 

where fij is the electrochemical potential of the species J and firef is the electrochemi- 
cal potential of the electrons in the reference electrode. When firef is changed to fi'^^j, 
with A/ire/ = /i^e/ — /^re/, the rcactiou (fT2|) will advance or retroced, modifying the con- 
centrations of J by negligible quantities but creating a macroscopic variation of (psoi, the 
electrostatic potential of the solution. If we call Acpsoi such a variation, elementary ther- 
modynamics gives: 

IJ,'j - fij = Afij = eCjA(j)soi (15) 
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The condition for the new thermodynamics equihbrium is given by Ujfi'j — nfi'^^j- = 0. 
Taking the difference between this equation and eq.( |T^ . one obtains: 

l^j^fij — nAflref = (16) 

j 

Using the expression of Afij as given in eq. lfTSil . eq. (fT6l) becomes: 

eA(l)soi ^ i^jCj - nAfiref = (17) 
j 

Exploiting the charge balance eq.( fT3|) . we finally get: 

- eA(j)sol = Afiref (18) 

B Numerical Methods 

In the article, we have briefly described the model used to calculate the electrostatic po- 
tential distribution (and thus aup, Osnfe and /?) for Azurin in the ECSTM setup. Here, 
technical details on these calculations will be given. We recall that we are using a contin- 
uum model to describe the system: the electrodes (tip and substrate) are assumed to be 
perfect conductors (the substrate is flat, the tip has a spherical curvature), the protein is 
a complex-shaped dielectric that cannot be penetrated by ions and the ionic solution is a 
dielectric that fills all the space left empty by the other components of the system. In par- 
ticular, the presence of an ionic atmosphere is taken into account by using the linearized 
Poisson-Boltzman equation eq.(2), which reduces to the Poisson one inside the protein. 
Thus, to calculate the electrostatic potential distribution in the system, we have to solve 
these equations with the boundary conditions given by the applied tip and substrate po- 
tentials. The numerical procedure used to perform such a task is the Finite Difference 
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(FD) method. Basically, the differential Poisson-Boltzmann equation is discretized on a 
regular grid, using standard expressions for the Laplacian (correct to the third order in 
the grid spacing, for the present article). Such a discretization translates the differential 
Poisson Boltzmann equation in a linear problem of the kind Ax = b, where the known 
term b is given by the electrostatic potential at the boundary of the considered system, 
the unknown x is the potential in the interior of the system and the coefficient matrix A 
depends on the dielectric constant e and the inverse screening length k. The boundary 
condition between the protein and the solution can be directly taken into account in the 
resolution of the problem by defining a position-dependent dielectric constant and inverse 
screening length in such a way that e(r) = Cpro and k{r) = if f G protein, and e(r) = €soi 
k{r) = ksoi if r ^ protein. We have defined the protein region as the interior part of the 
solvent-excluded surface of the protein. In particular, in all the simulations we have used 
^pro = 4, a typical value for electrostatic calculations in protein media. The step change of 
properties at the protein boundary is incompatible with the FD approach, where quantities 
should not vary too much on the grid spacing. For this reason, we smeared e(r) and k{r) 
at the boundaries by interpolating between the protein and the solution values of e and k. 
The smearing is obtained via the arctan function. 

To improve the description of the solution regions close to the electrodes (i.e., the 
Helmholtz layers), for the region outside the protein we have assumed that: 

e(f) = 6.0 k{f) = for d < 2.8A 
e(r) = 30.0 k(r} = k^oi for 2.8A <d< 4.8A 
e(f) = 78.39 k{r) = k^oi for 2.8A <d< 4.8A (19) 

where d is the distance of the point r from the surface of the closest electrode. The 
distances and the dielectric constants are typical values used to reproduce the experimental 
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capacitance of double layers. They are reminiscences, in the framework of our very simple 
model, of the properties of the inner {d ~ 2.8 A) and outer {d ~ 4.8 A) Helmholtz planes. 
Again, at the boundaries between the various d regions, e and k have been smoothly 
interpolated. 

At this point, we would like to discuss which boundary conditions we used for the 
electrostatic potential. We have discretized a rectangular box having the substrate and 
the tip surfaces as basis, and containing the protein at the center. On tip and substrate 
surfaces, we imposed a value of the ipsub/up potential given by (psub/up = Esub/up - E^ub/up^ 
taking (^soi = in eq.(2). On the lateral faces we imposed the potential distribution proper 
for the system without the protein. We enlarged the box until the potential results in the 
protein region did not depend on the box size anymore. 

We remark that the potential distribution for any combination of (psub and (pup values 
can be obtained with just two calculations: one with (psub = 1 and (pup = and the other 
with (psui, = and (pup = 0. Since the electrostatic problem that we solved is linear, 
the potential distribution for any other value of (psub and (pup can be obtained by linearly 
combining the results of these calculations. 

Finally, we give some details on the numerical resolution of the non-linear PB equation, 
necessary to obtain the results presented in Sect J4.3l An iterative procedure starting from 
the linearized solution were employed. At each step, the non-linear terms of the PB 
equation were Taylor-expanded to the first order in the potential, using the potential 
obtained in the previous step as the center of the expansion. The so-obtained linear 
equation is numerically solved and a new iteration is started. The procedure was stopped 
when the potential values changed less than a given relative threshold ( 10^^). To test the 
correctness of the results, we checked that the potential profiles presented in Ref.27 were 
reproduced by our code. 

We conclude this section by mentioning that the model described here has been im- 
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plemented in an home-made FORTRAN90 code, and that the solution of the numerical 
FD problem has been obtained by standard numerical techniques (bi-conjugate gradient 
algorithm). Typical parameters used in the calculations are: a grid space of 0.65 A, K 
300 grid points in the x and y directions (parallel to the electrode surface) and 60-80 grid 
points in the z direction (perpendicular to the electrode surface). 
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